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1.  Introduction 

Nanostructured  films  axe  widely  used  in  optics,  microelectronics,  laser  engineering  and 
many  other  fields  of  science  and  technology  [1],  The  materials  must  possess  certain 
properties,  for  instance,  mechanical  or  electrophysical  depending  on  the  application.  The 
properties,  as  a rule,  are  determined  by  the  composition  of  the  growing  films  and  also  by 
the  size  of  the  crystalline  grains  forming  the  films  [2].  Both  systems  can  be  studied  in 
a similar  way.  Film  systems  are  examined  here  only  as  an  example  since  now  they  are 
widely  used. 

Irrespective  of  the  preparation  techniques,  the  films  in  the  process  of  their  growth  pass 
a series  of  stages  typical  for  a first-order  phase  transition  [1,  3-9].  These  stages  are  the 
nucleation  of  nano-clusters  of  a new  phase  [1,  10-17],  their  growth  [10],  OstwaJd  ripening 
[1,  3-8,  17],  coagulation  [11]  and  formation  of  a continuous  film.  When  the  nucleation 
process  occurs  under  considerable  supersaturation,  the  density  of  islands  of  new  phase  and 
the  degree  of  substrate  filling  are  high  and  the  stages  following  nucleation  may  not  take 
place  [1].  In  this  case,  both  composition  and  structure  are  determined  by  the  nucleation 
processes  which  have  a complicated  non-linear  nature  [12-16].  To  control  the  process  at 
this  stage  is  not  possible.  Special  methods  such  as  nucleation  process  control  are  to  be 
used.  As  shown  in  Refs  [3-9],  the  growth  parameters  can  be  controlled  only  at  the  Ostwald 
ripening  stage,  where  nucleation  of  new  islands  does  not  take  place  and  a continuous  film 
has  still  not  formed.  The  possibility  of  composition  and  structure  control  is  based  on 
the  fact  that  at  this  stage  an  ensemble  of  islands  of  new  phase  begins  to  interact  with 
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generalized  diffusion  and  thermal  fields  formed  on  the  substrate’s  surface,  and  fluxes  of 
atoms  incident  on  it. 

The  aim  of  this  work  is  to  find  out  the  dependence  of  the  distribution  function  of  nano- 
islands on  their  composition  and,  accordingly,  the  distribution  functions  with  respect  to 
one  or  another  property  (such  as,  conductivity,  mobility  and  so  on)  dependent  on  the 
initial  composition.  It  was  shown  previously  that  there  is  a great  difference  between  evo- 
lution of  nano-islands  grown  from  substances  forming  a continuous  series  of  solid  solutions 
and  stoichiometric  substances  formed  as  a result  of  chemical  reactions.  This  leads  to  a 
significant  difference  in  the  growth  processes  of  films  both  at  the  nucleation  stage  and  the 
Ostwald  ripening  stage  [3-9].  For  this  reason,  these  cases  are  analyzed  separately. 


2.  Evolution  of  the  phase  composition  of  nano-islands  at  the 
Ostwald  ripening  stage. 


2.1.  EVOLUTION  OF  PHASE  COMPOSITION  IN  SOLID  SOLUTION 


Let  the  equilibrium  state  diagram  of  a system  for  a growing  film  be  the  state  diagram 
of  a binary  system  with  continuous  solubility  of  components  in  a solid  phase  (see  Fig  1). 
This  means  that  both  nucleation  and  Ostwald  ripening  processes  are  possible  only  when 
the  system  is  supercooled  to  a temperature  between  the  liquids  line  and  the  Ta  line  (see 
Fig  1).  When  the  system  is  cooled  in  such  a way  that  the  temperature  rapidly  arrives  at 
values  lying  between  the  Ta  — line  and  the  solidus  curve,  then  the  system  readily  solidifies 
without  separation  into  phases  and  the  composition  of  the  solid  phase  corresponds  to  that 
of  the  liquid  phase  [18].  These  processes  are  not  discussed  in  this  paper. 

If  a nano-island  ensemble  on  the  substrate  surface  is  at  the  Ostwald  ripening  stage  and 
is  supercooled  to  a temperature  in  the  interval  between  Ta  and  the  liquidus  line,  then  the 
composition  of  each  island  is  related  with  its  size,  R,  by  the  following  equation  [6,  7]: 


Csr  = 


kBTR 


+ Clo  + 


, f2<rwCm 
f V kgTR 


(1) 


where  <j  is  the  surface  tension  on  the  island-vapor  (solution)  interface;  kg  is  the  Boltz- 
mann constant  ; T is  the  absolute  temperature;  Cxo  is  the  equilibrium  concentration  of 
a component,  for  instance,  A on  the  liquidus  line  (Fig.l),  i.e.  in  a liquid  phase;  w is  the 
volume  per  atom  in  an  island;  rj)  is  the  functional  relationship  for  the  solidus-liquidus  line, 
i.e.  Cso  — Cxo  = V’(C'lo),  describing  the  diagram  of  state  (Fig.l);  Cs  — is  the  equilibrium 
concentration  of  a component  of  a solid  solution,  for  instance,  A in  the  solid  phase. 
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Figure  1:  Diagram  of  states  for  a binary  system  with  continuous  solubility  of  components  in  a solid  phase. 
Here  Cso  and  Clo  are  equilibrium  concentrations  of  a component  in  solid  and  liquid  phases,  respectively. 
Csr  is  the  composition  of  islands  of  radius  R\  Ta  is  the  temperature  below  which  an  alloy  solidifies 
without  changes  in  composition;  T(0)  and  Ci(0)  are  the  initial  temperature  and  concentration  in  melt 
and  solid  phase,  respectively;  a is  the  slope  angle  of  the  liquidus  line  for  C -4  0;  /3  is  the  slope  angle  of 
the  solidus  line  for  C -4  0. 


Let  us  consider  the  ranges  of  values  of  the  ip  function  near  small  values  of  Clo  concentra- 
tion of  the  component  A , i.e.  near  the  origin  of  coordinates.  In  this  case  the  ip  function 
can  be  approximated  by  two  straight  lines.  In  this  case 


Csr 


IffwCLO 

kBTR 


+ 


2owClo 

l^TR 


(2) 


Here  7 = tg  a/tg  /3  — 1;  tg  a is  the  angle  coefficient  of  the  straight  fine  binding  liquidus  line 
and  tg  P is  the  angle  coefficient  of  the  solidus  line;  (1  + 7)  is  the  distribution  coefficient 
of  liquid  and  solid  phases  within  the  equilibrium  one,  i.e.  1 + 7 = Cso/Clo . Note, 
that  if  7 < 0,  concentration  of  the  component  in  the  liquid  is  higher  than  in  the  solid 
(i.e.  an  easily  fusible  component  is  analyzed). In  the  case  when  (1  +7)  > 1,  i.e.  7 > 0, 
concentration  of  this  component  is  higher  in  the  solid  phase  than  in  the  liquid.  The 
expression  for  the  island  radius  R could  be  obtained  from  (2): 

2<7wCL0(1+7) 

kBT[cSR-cL0(i  + 'r)Y 

Let  us  define  Csr  — Clo(1  — 7)  as  fjj.  It  is  seen  from  equation  (3)  that  there  is  an 
unambiguous  relation  between  the  island  radius,  R,  and  composition.  Thus,  if  Csr  — > 
Cxo(l  + 7),  he.  £r  — > 0,  then  fi-4oo.  If  -4  00,  then  i?  — 4 0.  It  was  shown  previously 
in  Refs  [6-8],  that  in  order  to  find  out  the  size  distribution  function  of  islands,  f(R,t)  at 
the  Ostwald  ripening  stage,  it  is  necessary  to  solve  the  set  of  equations  that  involves  the 
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continuity  equation  for  the  distribution  function,  the  matter-and-heat  balance  equation, 
the  equation  describing  relation  between  the  flux  of  matter  arriving  at  an  island  and  the 
flux  of  heat  released  during  the  growth  of  the  island;  the  coupling  equation  of  the  melting 
point  for  an  island  with  radius,  R , with  the  concentration  of  matter  which  is  in  equilibrium 
with  it  and  the  equation  expressing  the  dependence  of  the  growth  rate  of  an  island  on  its 
radius. 

Since  there  is  an  unambiguous  correlation  between  R and  £r,  let  us  consider  the  distri- 
bution function  f(R,  t)  - the  distribution  function  of  island  versus  inverse  composition 
values,  i.e.  ip  In  accordance  with  (3),  for  t -t  oo,  and  R -t  oo,  (R  ->  0,  it  is  more 

convenient  to  use  the  reciprocals  of  the  compositions,  i.e.  1/£r  — » oo  at  t — t oo.  So,  we 
pass  from  f(R,t)  to  tp(l/£R,t),  i.e.  f(R,t)dR  = ip(l/£R,t)d(l/(R).  In  this  case  the  set  of 
equations  describing  the  size  distribution  function  f(R,t),  derived  in  Refs  [6,  7]  for  the 
function  <p(l /£R,t),  could  be  written  as: 


dt 


+ ms  hs'*1 


jO/gg) 

dt 


= 0 


QD(t)  = N0AC(t)  + * J <p(l/tR,t)(l/&W/tR) 
0 


QT(t)  = XCPPAT(t)  + «*/ ¥>(!/&,  *)(!/&)<*(  1/60 

0 


L Jd,r  = Jt,r 

Tr  = p(Cl,  R) 

<1(1/60  _ Ap  q(i/^fl) 

dt  ar  (l/ur1  Rc 


(4) 

(5) 

(6) 

(7) 

(8) 
(9) 


Here  Eqn.  (4)  is  the  continuity  equation  of  the  distribution  function  of  islands  versus 
island  composition;  Eqn.  (5)  is  the  equation  of  the  matter  conservation  law  expressed 
via  the  distribution  function  in  terms  of  composition;  Eqn.  (6)  is  the  equation  of  heat 
conservation  law  expressed  through  the  distribution  function  in  terms  of  composition; 
Eqn.  (7)  is  the  equation  coupling  heat  and  matter  fluxes  on  the  surface  of  an  island 
of  composition  Csr',  Eqn.  (8)  is  the  equation  coupling  the  melting  point  of  an  island  of 
radius,  R,  and  its  composition;  Eqn.  (9)  is  the  equation  describing  the  change  of  an  island 
composition  in  time;  Qo(t)  and  Qr(t)  are  the  power  carried  by  the  matter  flux  arriving 
at  the  substrate  and  the  power  in  the  heat  flow,  respectively;  A C{t)  and  A T(t)  are  the 
supersaturation  and  supercooling  on  a substrate  or  in  the  melt,  respectively  [6,  7];  Cp 
and  p are  the  heat  capacity  and  density  of  a melt  or  substrate,  respectively  depending  on 
where  the  ripening  process  takes  place;  A is  the  thickness  of  the  melt  or  substrate  layer, 
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where  the  Ostwald  ripening  process  takes  place,  A is  of  the  order  of  average  island  size, 
i.e.  A ~ R ( more  information  about  A value  and  its  calculation  can  be  found  in  [1]  ); 
Ap  and  p are  the  coefficients  dependent  on  the  mechanism  of  island  growth  [1,  6,  7];  the 
coefficient  p can  take  the  values  of  p = 2,  3,4  (in  the  case  of  evolution  of  faceted  islands 
p — 5);  Rc  is  the  critical  radius  of  islands;  Xcr  = K*  _ k/L]  L is 

the  latent  heat  of  crystallization  per  atom  of  the  precipitated  phase;  No  is  the  number  of 
adsorption  sites  on  the  substrate,  No  ~ 1/a2  (or  the  total  number  of  atoms  in  solution, 
Na  depends  on  the  kind  of  heat/mass  transfer  mechanism  realized  in  a system  [6,  7];  a 
is  the  lattice  constant  of  the  substrate;  %(#)  = 2(2  — 3 cos  8 + cos3#)-1;  6 is  the  wetting 
angle;  a =•  2awCu>{\  + 7).  Note,  that  in  the  general  case,  as  shown  in  Refs  [6,  7], 
the  values  of  Qo{t)  and  Qr(t)  can  be  represented  in  the  form  Qn{t)  — Qod  + 

Qp(t)  = Qqt  + gortn,  where  Qod  and  Qot  are  the  total  amount  of  matter  and  heat  on 
a substrate  at  the  beginning  of  the  Ostwald  ripening,  respectively;  goo  and  gQp  are  the 
strengths  of  sources  (sinks)  of  matter  and  heat;  n is  an  arbitrary  number,  n > 0 at  the 
Ostwald  ripening  stage,  i.e.  when  t — » oc,  Qd( f)  — > gootn  and  Qp(t)  — »■  poT<"- 
Equations  (4)-(9)  were  derived  by  substitution  of  variables  from  Eqn.  (3)  into  the  equa- 
tions describing  evolution  of  the  size  distribution  function  [6,  7].  Since  Rc(t)  — where 


— ct-Cr 

Cro(  1+7)  > 


then 


a/tR  _ (C{t)-CL o)(l+7) 
Rc(t)  CR-CL0(l+y) 

As  <7*0(1  + 7)  = Cso,  then 

q _ (C(t)  - CL0)CS0 

(rRS)  (CR  - C.so)CLO 


(10) 


(11) 


In  accordance  with  (11)  we  can  rewrite  (9),  as 


d(i/60  \{C{t)-cLo){  1 + 7) 

dt  (1  /^)p_1  . CR  — Cso 


(12) 


where  A'p  = Apj a. 

Finally  the  set  of  Eqns  (4)-(8)  and  (12)  completely  describes  the  evolution  of  phase 
composition  at  the  stage  of  the  Ostwald  ripening  of  a substance  with  the  state  diagram 
of  a binary  system  with  continuous  solubility  of  components  in  a solid  phase.  Let  us 
introduce  the  following  variables:  U = ’ x = Nit)’  T = !n a=P-  where  A0  and 

A (t)  are  the  supersaturation  at  the  initial  stage  of  ripening  process  n at  time  t.  Equations 
(4)-(9)  written  with  these  variables  are  completely  identical  to  the  equations  describing 
the  process  of  non-isothermal  Ostwald  ripening  [6-8].  Their  solutions  in  variables  u and 
r will  be  the  same.  The  solution  of  the  equations  are  given  in  [6-8].  By  passing  from 
variables  u and  r to  variables  1/£r  and  t and  also  to  the  function  ip(l/£R,t),  we  shall 
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obtain  the  relationship  of  island  distribution  with  respect  to  composition  in  time.  This 
solution  is  essentially  dependent  on  a technique  of  deposition  of  material  on  a substrate 
as  well  as  on  the  intensity  of  cooling  (or  heating)  of  the  substrate.  It  is  well  known, 
that  the  Ostwald  ripening  stage  starts  only  in  the  presence  of  weak  sources  of  matter  at 
0 < n < 3/p.  In  this  case,  as  analysis  shows  [6,  7]. 


¥>(!/&,*) 


2auCL0{l+'r)N{t) 

kBTRc(t)  A ’ 


(13) 


Here,  Pp(U ) is  the  probability  density  for  U to  range  between  U and  ( U + dU)\  the 
index  p indicates  the  particular  mechanism  of  heat/mass  transfer  realized  in  the  system; 
N(t)  is  the  density  of  islands  on  a substrate  at  the  time  t.  the  analytical  expressions 
of  this  value  are  given  in  [1,  3-8];  Rc(t ) is  the  critical  size  of  islands  at  the  time  t.  In 
Eqn.  (13)  instead  of  the  critical  radius  Rc{t)  one  can  use  the  average  radius,  since  they 
are  related  to  each  other  by  equation  R(t)  = Rc(t)Cpn,  (see  Refs  [1,  3-5]).  Coefficient, 
Cpn,  depend  on  the  mechanism  of  heat/mass  transfer  and  the  power  carried  by  sources 
of  matter  and  heat.  Numerical  values  of  Cpn  can  be  found  in  Refs  [1,  3-5].  The  value 
= AWis  supersaturation. 

The  value  of  supersaturation  in  accordance  with  [7]  changes  in  time  as: 


A(t)  = 


‘2aujCL0 

hBTCpn 


(14) 


Here,  Ro  is  the  average  size  of  islands  at  the  initial  stage  of  the  Ostwald  ripening.  Coef- 
ficient A®  is  the  kinetic  coefficient  associated  with  the  coefficient,  Ap,  derived  above  and 
it  depends  on  the  heat/mass  transfer  mechanism  realized  in  the  system.  Values  of  these 
coefficients  for  all  possible  mechanism  of  heat /mass  transfer  are  given  in  [6,  7].  In  partic- 
ular, when  the  islands  are  in  the  form  of  spherical  segments  and  their  growth  occurs,  as 
a rule,  due  to  heat /mass  transfer  on  the  substrate  surface  with  A,  R\  p = 4 
,0  27  Da(7w2N0CsboKs£3b(f>i(9) 

p ~ 32{DaN0CLMLHn(H/R1)  + Ks£,bkBT2ln{Xa/R1)y 


where  As  is  the  free  path  length  of  adatoms;  Da  is  the  Brownian  diffusion  coefficient 
of  adatoms;  <7  is  the  surface  tension  at  the  old  phase-island  interface  (solution,  melt, 
vapor);  </i (0)  is  the  geometrical  factor;  Ks  is  the  thermal  conductivity  of  the  substrate, 
isTjfJm-1-1];  t,b  is  the  thickness  of  a substrate  layer  where  the  ripening  process  takes  place 
[1,  7];  L is  the  latent  heat  of  crystallization  per  atom  of  a precipitated  phase;  CBsb o is  the 
equilibrium  concentration  of  adatoms  on  a substrate;  To  is  the  equilibrium  temperature 
of  crystallization  (of  solution,  vapor  or  melt);  H is  the  average  distance  between  islands 
on  a substrate  [7];  Ri  = R ■ sin$;  0 is  the  wetting  angle. 

Analytical  expressions  of  Pp(U)  functions  for  all  values  of  p are  given  in  [1,  3-8],  that  is 
why  in  Fig  2 only  the  plots  of  Pp(U ) values  for  p = 2, 3 and  4 for  n = 0 are  presented. 
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Figure  2:  Pz{U)^  P$(U)  and  P*(U)  functions  for  nano-islands  with  the  spherical  shape. 


Pp(U)  functions  can  be  experimentally  found  using  the  technique  described  in  [1],  i.e. 
by  expressing  U via  the  average  radius  of  islands  U — R • Cpn/ R.  This  way  is  more 
convenient  than  to  express  U via  concentrations  U = (C(t)  — Csq)/(Cr  — Cs o),  which 
are  hardly  measurable  values  in  the  experiment.  The  distribution  function  in  terms  of 
composition  (13)  can  be  actually  expressed  through  values  independent  of  composition 
(except  the  Cs o value)-this  is  an  advantage  of  our  approach.  Equilibrium  concentrations 
Clq  and  Ca o as  well  as  equilibrium  temperature  T0  are  determined  from  the  solution  of 
Eqns  (5)-(8)  for  t -4  oo.  For  t o o the  equations  transform  into  algebraic  equations. 
The  equations  can  be  easily  solved  and,  as  a result,  the  values  of  Clo,  Cso  and  To  can  be 
obtained  (for  more  details  see  [6,  7]). 

The  average  composition  of  nano-islands  is  calculated  from  the  equation 


m = 


oo 

J 


(nit) 


( 1 ^ 
V6?  j 


(15) 


Taking  into  account  Eqn.  (3)  and  expressing  p{l/^R,t)d(l/^R)  via  f(R,t)dR,  as  it  is 
done  before,  we  shall  get 


m 


2(7a>Cto(l  + 7) 


kBT 


my1  ■ 


Finally  we  shall  obtain: 


l%  + A°rt}-'lr 


(16) 


(17) 


Comparing  Eqn.  (17)  with  (14),  A (t)  = £(f)/(l  +7)  is  readily  obtained.  This  result  is 
associated  with  the  fact  that  compositions  of  solid  and  liquid  phases  are  related  to  each 
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Figure  3:  Dependence  of  composition  of  nano-islands  on  size  (a)  and  composition  of  nano-island  on  time 
(b);  time  dependence  of  the  mean  composition  of  nano-islands  (c). 

other  by  the  state  diagram  presented  in  Fig.l. 

The  changes  in  composition  in  an  island  can  be  obtained  from  the  following  relationship 

CR  = CS0  + ^1,  (18) 

where  A (f)  is  determined  from  Eqn.  (14),  whereas  U = R/Rc. 

The  dependence  of  island  compositions  on  their  sizes  are  shown  in  Figures  3a,b,c  (for 
composition  of  the  particular  island  dependent  on  its  growth  time  and  dependence  of  the 
average  composition  of  islands  on  time). 

Let  us  do  some  estimates.  As  follows  from  Eqn.  (18).  Cr  = Cso  + 2£bTR  • us  assume, 
a to  be  ~ 10-1  -r  10 J/m2,  and  kT  ~ 10~19J,  we  shall  get  for:  Cr  ~ Cso  (l  + — ) ■ For 
this  reason,  the  main  change  in  composition  of  islands  takes  place  for  R < 10_8m.  In  this 
case,  the  change  of  composition  will  be  within  a few  tens  of  atomic  percents  and  this  fact 
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is  of  particular  importance  for  the  measurement  electrophysical  properties. 

2.2.  EVOLUTION  OF  THE  COMPOSITION  OF  A NANO-ISLANDS  OF  STOICHIO- 
METRIC COMPOSITION. 

When  islands  are  of  strictly  stoichiometric  composition,  as  it  is  shown  in  Refs  [1,  3-5,  8], 
the  regions  of  phase  coexistence  and  accordingly,  the  composition  of  island  films  are 
determined  from  the  solution  of  the  system  of  equations  derived  in  [1,  5].  Detailed  analysis 
of  the  set  of  equations  is  performed  in  Ref.  [1,  3-5,  8]  and  some  particular  solutions  for 
specific  systems  are  given.  Now  we  shall  pass  to  the  description  of  phase  composition 
evolution  at  the  nucleation  stage. 


3.  Evolution  of  phase  composition  in  new  phase  islands  at  the 
nucleation  stage 

3.1.  EVOLUTION  OF  PHASE  COMPOSITION  IN  SOLID  SOLUTION  ISLANDS. 

At  the  nucleation  stage  there  is  no  unambiguous  relation  between  size  and  composition 
as  compared  to  the  Ostwald  ripening  stage.  Note,  that  new  phase  nucleation  processes  in 
such  systems  were  investigated  by  many  authors  [12,  14,  15,  20,  21].  It  turned  out  that 
fluctuations  of  components  involved  in  a nucleus  happen  independently  of  each  other. 
The  control  of  the  composition  is  impossible  at  this  stage. 

3.2.  EVOLUTION  OF  PHASE  COMPOSITION  IN  ISLANDS  OF  STOICHIOMETRIC 
COMPOSITION. 

When  islands  of  new  phase  are  stoichiometric  compounds,  the  nucleation  process  can 
occur  in  several  stages  (see  Ref.  [20]).  It  is  shown  that  the  nucleation  process  is  essentially 
dependent  on  whether  or  not  chemical  reactions  in  multicomponent  adsorbed  gas  occur 
simultaneously  with  the  formation  of  the  end  reaction  product.  We  have  shown  in  [20] 
that  for  specific  conditions  it  is  possible  to  observe  self-oscillations  of  the  number  of  nuclei 
and  self-organization  phenomena  are  possible. 
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4.  Evolution  of  the  composition-dependent  properties  of 
growing  nanostructured  films 

It  is  well  know  that  many  properties  of  solids  are  essentially  dependent  on  their  chemical 
and  phase  composition.  The  properties  can  be  the  concentration  of  charge  carriers  in 
semiconductors,  their  mobility,  conductivity  of  metals  and  alloys,  absorption  and  electro- 
magnetic wave  scattering,  strength  properties  [2, 18],  etc.  Without  claiming  any  consistent 
investigation  of  the  properties,  we  show  here  merely  the  main  principle  for  calculation  of 
their  evolution  in  the  process  of  nanostructured  film  growth. 

4.1.  EVOLUTION  OF  PROPERTIES  IN  SOLID  SOLUTION  NANO-ISLANDS. 

Let  F(C,  t)  be  the  function  describing  the  relation  between  certain  properties  of  the  solid 
(in  this  case,  an  island)  on  composition.  Let  us  find  the  distribution  function  of  properties 
of  different  composition.  To  do  this  we  shall  pass  from  the  function  ^>(1  /£r),  expressing 
distribution  function  of  island  in  terms  of  composition  to  the  distribution  function  with 
respect  to  properties,  & (F,  t): 

&(F,t)d(F)  = <pQ-,  (19) 

In  the  general  case,  the  average  value  of  properties  in  nano-islands  changes  in  time  as 

oo 

F(t)  = / F(t)Q(F(t),t)d(F)  . (20) 

o 

Note,  that  Eqns  (4)-(9)  describe  not  only  the  evolution  of  composition  and,  accordingly, 
properties  in  nano-island  films  but  the  meaning  of  these  values  in  continuous  films  as  well. 
For  this  purpose  it  is  necessary  to  estimate  the  time  of  coalescence  of  islands  into  a solid 
layer  using  the  equations  given  in  Ref.  [3]  and  the  distribution  functions  of  islands  in  terms 
of  composition  in  time.  In  this  case,  it  is  necessary  to  take  into  account  that  Eqn.  (4)-(9) 
and,  accordingly,  ( 13)— ( 18)  are  valid,  strictly  speaking,  only  for  total  coverage  close  to 
zero.  In  the  case  when  the  total  coverage  essentially  exceeds  zero,  the  collision  integral  is 
to  be  introduced  into  Eqns  (4)-(9),  [22],  However,  as  it  follows  from  calculations  [22],  the 
function  values  are  only  slightly  changed.  Therefore,  Eqns  (13)— (18)  provide  an  adequate 
approximation. 

When  some  property  of  islands  depends  both  on  their  size  and  composition,  i.e.  0 = 
0(C,  R.t),  it  is  necessary  to  pass  from  variable  R in  the  terms  for  0(U,  R,  t)  to  variable 
Csr  by  means  of  Eqn.  (3).  In  a consequence,  we  again  get  the  function  0(C,  t).  In 
the  case,  when  the  property,  in  question,  is  only  the  distribution  function,  i.e.  0(R), 
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conditions  for  its  evolution  can  be  expressed  through  the  size  distribution  function  f{R,t) 
[!]• 

4.2.  EVOLUTION  OF  PROPERTIES  IN  NANO-ISLANDS  OF  STOICHIOMETRIC 
COMPOSITION. 


The  nano-islands  of  stoichiometric  composition  are  homogeneous.  Therefore,  the  property 
function  depends  on  the  islands’  size  only.  The  properties  of  a nano-film  are  determined 
by  the  number  of  islands  of  one  or  another  phase  and  by  the  size  distribution  of  islands 
colliding  with  each  other.  The  technique  for  the  phase  coexistence  regions  determination 
is  described  in  detail  in  Ref.  [1,  3-5].  The  diagram  of  phase  coexistence  regions  for  a three- 
component  two-phase  system  (where  one  of  the  components  is  simultaneously  included 
in  both  phases)  is  constructed  therein.  If  there  is  an  unambiguous  relation  between  the 
distribution  function  of  property  values  F for  s-phases,  Qs(F,t),  and  size  of  particles, 
then  the  distribution  function  for  n < 3/p  can  be  derived  in  a similar  way  as  in  (13): 


Qs(F,t) 


N‘(t)  dR 
Ri(t)  dF‘ 


p;(u), 


(21) 


where  N‘(t ),  P‘(U),  R‘(t)  are  the  functions  of  the  number  of  islands  Ns(t),  probability 
density  P‘{U)  and  critical  radius  of  islands  of  .s-phasc  [1,  3-5];  dFa/dR  is  the  derivative  of 
relation  between  the  properties  and  radius  of  an  island  of  s-phase;  p = 2, 3,4  in  accordance 
with  the  island  growth  mechanism.  Evolution  of  the  average  value  of  a property  can  be 
described  by  an  equation  similar  to  Eqns  (15). 

Consequently,  all  conclusions  on  influence  of  heating  or  cooling  of  a substrate  and  the 
power  of  component  sources  presented  in  Refs  [1,  3-5]  for  the  size  distribution  function 
f’(R,t)  of  islands  of  phase  s over  a substrate  can  be  extended  for  the  function  Qs(F,t). 
When  islands  of  different  phases  (chemical  compounds)  have  regions  of  mutual  solubility, 
i.e.  they  can  be  solid  solutions,  evolution  of  their  properties  can  be  described  in  a similar 
way  as  in  section  2. 

On  the  contrary,  when  properties  depend  on  sizes  of  islands  and  the  function 
Q(F,t)  is  known,  the  relative  number  of  s-phase  can  be  expressed  as  Ja  = 
f fs{R,t)R3dR;  (here  V’n  is  the  volume  of  molecules  of  phase  s:  N0  is  the 

number  of  adsorption  sites  on  the  substrate  surface  N0  ~ 1/a2;  a is  the  lattice  constant 
of  a substrate)  can  be  expressed  as  follows: 


Js  = 


1 


k(6)  I e°(F,t)Z3(F)d(F), 

0 


(22) 


taking  into  account  the  condition  f”(R,t)dR  = Q’(F,  t)dFs.  R = Z(F)  is  the  inverse 
function  from  the  function  F = Z(R).  Thus,  the  coexistence  region  of  different  compo- 
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Figure  4:  Values  diagram  of  properties  in  the  space  of  components  Qi  and  temperature  ( in  projection  on 
the  Q i,  Q 2 ) plane.  Boundaries  of  the  values  of  properties  in  a system  are  indicated  with  the  continuous 
line;  displacement  of  the  boundaries  due  to  temperature  change  is  indicated  with  dashed  and  chain-dotted 
line.  Q°  and  Q°  are  certain  initial  quantities  of  components;  I and  II  are  the  regions  of  different  phases; 
F1  and  F1 1 are  the  values  of  a property  in  islands  of  different  chemical  composition.  The  trajectory  of 
a system  in  the  presence  of  sources  of  components  is  indicated  with  arrow. 

sitions  can  be  determined  by  known  values  of  F.  The  diagram  of  values  of  the  property 
function  for  a two-phase  tree-component  system  where  one  of  the  components  is  simul- 
taneously included  into  both  phases,  is  given  in  Fig. 4.  The  diagram  is  calculated  in 
accordance  with  the  technique  described  in  [1,  3-5].  Fig.  4 shows  that  in  the  case  when 
islands  form  a continuous  film,  the  structure  with  a number  of  properties  determined  only 
by  its  dependence  on  the  size  of  islands  exists  in  regions  I and  II.  Within  the  region 
/ + II,  islands  of  mixed  composition,  i.e.  with  a different  value  of  properties  F1  and  F11 
will  grow.  In  the  process  of  growth,  the  islands  give  rise  to  non-homogeneous  properties  of 
thin  films.  The  regions  of  F values  can  be  varied  in  the  presence  of  sources  of  components 
(or  through  cooling  heating  of  the  substrate).  For  this  purpose,  it  is  necessary  to  choose 
a certain  value  of  the  decay  degree  of  sources,  n,  in  accordance  with  the  island  growth 
mechanism  [1,  3-5], 

Now,  we  pass  to  the  nucleation  of  new  phase  islands  in  multicomponent  systems  during 
chemical  reactions.  It  is  shown  in  Ref.  [21]  that  for  some  specific  relation  between  chemical 
components  and  islands  of  different  phases  complex  non-linear  interactions  occur.  This 
causes  periodic  change  of  the  number  of  islands  in  space  and  time.  As  a result  periodic 
changes  in  thin  films  properties  are  observed.  The  value  will  change  periodically  in  space 
and  time  as  well  (see  Fig.  5a  and  5b). 

In  the  case  of  the  formation  of  islands  of  several  phases  (chemical  compounds),  the  be- 
havior of  properties  is  more  complicated.  However,  when  self-organization  or  self-induced 
oscillations  occur  in  a system,  appropriate  properties  will  also  periodically  change  either 
in  space,  or  in  time. 
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5.  Discussion 


The  relationship  between  size  and  composition  of  new  phase  nano-islands  formed  in  the 
course  of  condensation  of  solid  solution  films  is  elucidated  in  the  work.  Analysis  shows 
that  the  most  favorable  ways  to  control  the  composition  of  the  nanostructured  films  are 
at  the  Ostwald  ripening  stage.  At  this  stage  there  is  an  unambiguous  relation  between 
the  size  distribution  function  of  islands  and  distribution  function  in  their  properties.  This 
allows  to  pass  from  equations  describing  evolution  of  the  size  distribution  function  to 
equations  of  evolution  of  distribution  function  of  properties  in  solid  solution  nano-islands. 
The  analysis  shows  that  considerable  changes  in  composition  occur  for  sizes  less  than 
10_8m.  In  this  case,  their  composition  can  be  changed  by  tens  of  atomic  percent  that 
corresponds  to  1012  4-  1013m~3.  For  values  of  island  radius  R > 10 ~sm  the  composition 
of  islands  is  independent  of  size  actually  regardless  of  substance. 

It  is  shown  that  the  coexistence  diagrams  can  be  constructed  for  stoichiometric  com- 
pounds, and  by  means  of  change  of  atomic  fluxes  and  controlling  the  degree  of  cooling  or 
heating  of  a substrate  one  can  pass  from  one  property  to  another. 

It  is  found  that  the  control  of  properties  at  the  nucleation  stage  is  restricted.  Only  in 
the  case  with  nucleation  of  stoichiometric  compounds  the  structures  with  given  properties 
may  be  obtained,  provided  that  self-induced  oscillations  and  self-organization  take  place. 
However,  in  this  case  it  is  necessary  either  to  calculate  beforehand  the  oscillation  period 
of  the  composition  [21],  or  to  determine  it  experimentally  and  to  stop  the  growth  when 
the  desired  composition  is  achieved. 

This  work  was  partially  supported  by  the  Russian  Foundation  for  Basic  Research  (code  99- 
03-32768  and  98-03-32971),  grant  ’’Integration”  No.  A0121,  NATO  Grant  SfP  #973252. 
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